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Further Constructions of Control-Lyapunov Functions 
and Stabilizing Feedbacks for Systems Satisfying 
the Jurdjevic-Quinn Conditions 
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Abstract 

For a broad class of nonlinear systems, we construct smooth control-Lyapunov functions whose 
derivatives along the trajectories of the systems can be made negative definite by smooth control laws 
that are arbitrarily small in norm. We assume our systems satisfy appropriate generalizations of the 
Jurdjevic-Quinn conditions. We also design state feedbacks of arbitrarily small norm that render our 
systems integral-input-to-state stable to actuator errors. 
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I. Introduction 

Lyapunov stability is of paramount importance in nonlinear control theory. In many important 
applications, it is very beneficial to have a continuously differentiable Lyapunov function whose 
derivative along the trajectories of the system can be made negative definite by an appropriate 
choice of feedback. Observe in particular that: 

• Recent advances in the stabilization of nonlinear delay systems (e.g., [7], [13], [21]) are based 
on knowledge of continuously differentiable Lyapunov functions. 

• Lyapunov functions are very efficient tools for robustness analysis. For example, many proofs 
of nonlinear disturbance-to-state L p stability properties rely on Lyapunov functions; see [6, 
Chapter 13] and [3], [10], [15]. Moreover, control-Lyapunov function (CLF) based control 
designs guarantee robustness to different types of deterministic [5] and stochastic disturbances, 
and to unmodeled dynamics [16], [17]. 

1 Corresponding Author: F. Mazenc. The first author was supported by the MERE Project. The second author was supported 
by NSF Grant 0424011. F Mazenc is with the Projet MERE INRIA-INRA, UMR Analyse des Systemes et Biometrie INRA, 
2, pi. Viala, 34060 Montpellier, France, mazenc@helios.ensam.inra.fr. M. Malisoff is with the Department of Mathematics, 
Louisiana State University, Baton Rouge, LA 70803-4918, malisoff@lsu.edu. 
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• When a CLF satisfying the small control property (as defined below) is available, the universal 
formula in [19] provides an explicit expression for a stabilizing feedback that is also an optimal 
control for a suitable optimization problem whose value function is the CLF; see [19]. 

• Backstepping and forwarding require Lyapunov functions of class C 1 for the subsystems [17]. 
The converse Lyapunov theorem (see [9]) ensures that, for any system that is globally asymp- 
totically stabilizable by C 1 feedback, a CLF exists. Unfortunately, for nonlinear control systems, 
determining explicit expressions for CLFs is in general difficult. Fortunately, for large classes of 
systems, one can determine functions whose derivatives along the trajectories can be rendered 
negative semz'-definite. If the systems satisfy the so-called weak Jurdjevic-Quinn conditions 
(defined below), which generalize those given in [8], then globally asymptotically stabilizing 
feedbacks can be constructed. However, in this case, explicit formulas for CLFs are not generally 
available. This motivates the following fundamental question: When the Jurdjevic-Quinn method 
applies, is it possible to design explicit CLFs? 

In [4], where this issue was addressed for the first time, a method was presented for designing 
explicit CLFs for affine homogeneous systems that satisfy the Jurdjevic-Quinn conditions. Our 
objective in the present note is to extend the main result of [4] by constructing CLFs for systems 
satisfying appropriate generalizations of the Jurdjevic-Quinn conditions, but not necessarily 
having the homogeneity property, including cases where the system may not be control-affine. 
Our work also complements [14] where strong Lyapunov functions are constructed for a large 
family of systems satisfying either an appropriate Lie algebraic condition or which can be shown 
to be stable using the LaSalle invariance principle. The main difference between the present work 
and [14] is that in [14], only systems without input are considered whereas here we consider 
systems with input. 

We end this introduction by recalling some basic facts on the Jurdjevic-Quinn method. We 
say (see for example [4] for the relevant definitions) that a nonlinear control-affine system 

x = f(x) +g(x)u , g(x) = (gi(x), . . . ,g m (x)) (1) 

satisfies the (weak) Jurdjevic-Quinn conditions provided there exists a function V : W n — > R 
satisfying the following three properties: V is positive definite and radially unbounded; for all 
x E M n , LfV(x) < 0; and there exists an integer / such that the set 

W(V) = [x e W 1 : Vfc e {1, . . . , m} and i e {0, . . . , /}, L f V(x) = L ad)[gk) V(x) = o} 

February 2, 2008 DRAFT 



3 



equals {0}. Here and in the sequel, we assume all functions are sufficiently smooth. If ([TJ 
satisfies the weak Jurdjevic-Quinn conditions, then it is globally asymptotically stabilized by 
any feedback u = —£(x)L g V(x) T where £ is any positive function of class C 1 . The proof of 
this result relies on the LaSalle Invariance Principle. 

The remainder of this paper is organized as follows. In Section HH we present our main result. 
Section |Inl is devoted to a discussion of our main result, Section |E1 to its proof, and Section Ivl 
to an illustrating example. Section |VI| constructs feedbacks for our systems that have arbitrarily 
small norm and that in addition achieve integral-input-to-state stability relative to actuator errors. 
Concluding remarks in Section IVIII end our work. 



II. Main result 

Recall (cf. [2]) that a C 1 positive definite function V(-) on M n is called a control-Lyapunov 
function (CLF) for a system x — Piix) + L P2(x) u with input u provided it is radially unbounded 
and satisfies L vl V(x) > => [x = or L tp2 V(x) 0]- We use V(x,u) to denote the 
derivative V(x, u) = L ipi V(x) + L ip2 V(x)u of V along trajectories of the system. We often 
suppress the arguments of V to simplify the notation. We say that a CLF V(-) for the system 
X = fxix) + L P2{x) u satisfies the small control property [19] provided for each e > 0, there 
exists 5(e) > such that if < \x\ < 5(e), then there exists u (possibly depending on x) sucn 
that \u\ < e and L ipi V(x) + L ip2 V(x)u < 0. 

We next provide our main CLF and stabilizing feedback constructions for the fully nonlinear 
system 

x = F(x, u) (2) 

where x G W 1 , u G M. m is the control, F(0, 0) = 0, and the function F is assumed to be C 1 . We 
further assume that u i— > F(x, u) is C 2 (i.e., the second order partial derivatives, with respect to 
the components of u, of each component of F are continuous), so the functions 

f(x):=F(x,0), g(x):=^{x,0) (3) 

are at least C 1 . Finally, we assume: 

Assumption 1: A smooth function V(x) that is radially unbounded and positive definite and 
such that 

LfV(x) < Vxel" (4) 
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is known. 

Assumption 2: A vector field G(x) such that if L g V(x) = and x ^ 0, then we either have 
L f L G V(x) < or LjV(x) < 0, is known. 
We are ready to state our main result. 

Theorem 1: Assume the data © satisfy Assumptions [UE1 Then one can determine a positive 
definite smooth function 5 : [0, oo) — > [0, oo) and a function Cl : [0, oo) — > [0, oo) such that 

r V(x) 

V\x) = V(x)+ / n(s)ds + 8(V(x))L G V(x) (5) 
Jo 

is a CLF for © that satisfies the small control property. In fact, for each real-valued C 1 positive 
function £(•), one can determine a function S(-), and a C 1 function £ : [0, oo) — > (0, oo) satisfying 
< f° r s > 0, such that © is a CLF for © satisfying the small control property 
whose derivative along the trajectories of in closed-loop with the feedback 

u = -t(V(x))L g V(x) T (6) 

is negative definite. 

III. Discussion of Theorem [Tj 

1. Assumptions [T] and El are similar to the assumptions of the main result of [4]. In particular, for 
the special case where F is control-affine, [4] provides an explicit expression for a vector field 
G(x) such that Assumption |2] holds whenever the so-called "weak Jurdjevic-Quinn conditions" 
(see the introduction) are satisfied. This vector field is not continuous at the origin but it turns 
out that there exists an integer N > 1 such that the vector field Gj^(x) = V(x) N G(x) is of class 
C°° for V satisfying our assumptions. The equality LjLg n V(x) = NV (x) N ~ l L fV (x)LgV (x) + 
V(x) N LfLaV(x) then implies that if G(x) satisfies Assumption El and if Assumption [T] also 
holds, then Gn(x) satisfies Assumption |2] as well. Consequently, one can take advantage of the 
formula in [4] to determine a C°° vector field for which Assumption El is satisfied. 

2. No restriction on the size of the function £(•) in © is imposed. Therefore, the family of 
feedbacks © contains elements that are arbitrarily small in (sup) norm. In fact, for any continuous 
positive function e : [0, oo) — > (0, oo), we can design our stabilizing feedback u so that it satisfies 
\u{x)\ < e{\x\) for all x G K". 



February 2, 2008 



DRAFT 



5 



3. An important class of dynamics covered by Theorem [T] is described by the so-called Euler- 
Lagrange equations 

d fdL . \ dL . .. 

JtW q - q) )-^ ( "'" ) = T <7) 

for the motion of mechanical systems, in which q represents the generalized configuration 
coordinates, L = K — P is the difference between the kinetic energy K and potential energy 
P, and r is the control [22]. In standard cases, K(q,q) = \q T M(q)q where the inertia matrix 
M(q) is C 1 and everywhere symmetric and positive definite. Then the generalized momenta |p 
are given by p = M(q)q, so in terms of the state x = (q,p), the equations © become [22] 

dH dH 

q= P) T = M- 1 (q)p J p=- ir (q,p) T + t, (8) 

op oq 

where H(q,p) = ^p T M~ 1 (q)p + P(q) is the total energy of the system. We make the following 
additional assumptions: (a) P(q) is positive definite and radially unbounded and (b) VP(q) ^ 
whenever q ^ 0. (These two assumptions are not too restrictive since one can often modify H 
and r to get a new system that satisfies these assumptions. Condition (a) can be weakened by 
assuming there is a constant c such that q i— > P(q) + c is radially unbounded and positive definite 
in which case we simply add c to the function V in what follows.) Then H is positive definite and 
radially unbounded, so V = H satisfies Assumption [T] The radial unboundedness follows from 
the continuity of the (positive) eigenvalues of the positive definite matrix M -1 (g) as functions of 
q [20, Appendix A4], which implies that each compact set S of q values admits a constant c$ > 
such that p T M~ 1 (q)p > cs\p\ 2 for all q E S and all p. In our general notation with x = (q,p), 
we get L f V(x) = and L g V{x) = H p (x) = p T M-\q). Choosing G{x) = [ VP(g) T ] T 
gives L G V(x) = H p (x)VP(q) T . Therefore, if L g V(x) = p T M- 1 (q) = and x ^ 0, then p = 
and therefore also LjLgV(x) = — V ' P(q)M~ l (q)V P(q) T and q ^ 0. Since M _1 is everywhere 
positive definite, Assumption |2] therefore reduces to our assumption (b) and therefore is satisfied 
as well. We study a special case of ® in Section |V] below, where we explicitly compute the 
corresponding CLF and stabilizing feedback. 

IV. Proof of Theorem [Tj 

A. Control Affine Case 

We fix a positive function £ : [0, oo) — > (0, oo), and functions V and G satisfying Assumptions 
HE We begin by proving Theorem [l] for the case where is control affine, i.e., of the form 
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dU). In this control affine case, the conclusions of our theorem will hold with f2 = and £ = £ . 
In Section IIV-B1 we will modify our constructions to handle the fully nonlinear system ©. 
First step. We exhibit a family of functions 5(-) for which the function 

U(x) := V(x) + 8{V{x))L G V{x) (9) 

is positive definite and radially unbounded. One can determine ot f (-) of class K,^ such that 

tti(kl) < < a 2 {\x\) and |L G V(x)| < a 3 (\x\) for all x G R n . It follows that 

Z7(ar) > ai(|z|) - 8(V(x))a 3 (\x\) > a 1 {a^ l {V{x))) - 5(V(x))a 3 (a^(V(x))) (10) 
for all x E W l . We can use standard results to find a C 1 function 5 : [0, oo) — > [0, oo) such that 

m < . ;ff ( :» v„> . (ii) 

With such a function £(■), the inequality U(x) > |ai(o;^ :L (y(s))) for all iGl" follows from 
(TTOb . Since V(x) is positive definite and radially unbounded and ls °f class /Coo, 

this implies that {7 (#) is positive definite and radially unbounded as well. In the next steps, we 
impose further restrictions on S. 

Second step. Along the trajectories x(t) of our system (HJ) in closed-loop with the feedback 
u = —£(V(x))L g V(x) T , the derivative U of U(x) from © reads 

U = [L f V(x) - aV(x))\L g V(x)\ 2 ] [1 + 5'(V(x))L G V(x)] 

+ 5(V(x))L f L G V(x) - £(V(x))6(V(x))L g L G V(x)L g V(x) T . (12) 

We restrict our attention to functions 5 such that 

5'(V(x))L G V(x) > -~ ViG M n . (13) 

Recalling © and (fT2l therefore gives the inequality 

U < \[L f V{x)-l{V{x))\L g V{x)\*] + 5(V(x))L f L G V(x) 
+£(V(x))6(V(x))\L g L G V(x)\\L g V(x)\. 
^From ©, we deduce that 



(14) 



U < l -[L f V{x)-l{V{x))\L g V{x)\ 2 ]+5{V{x))L } L G V{x) 

+£{V(x))6 2 {V(x))\L g L G V{x)\ 2 . (15) 
Third step. The remaining part of the proof relies extensively on the following: 
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Lemma 2: Assume that the system © satisfies Assumptions [T]|2j Then, there exist continuous 
positive definite functions Y and iV satisfying the following: If |L ff V(a;)| < r(|x|), then either 
L f V(x) < -N(\x\) or L f L G V(x) < -N(\x\). 

Proof: We first show that the continuous function 

S(x) = mm{0,L f L G V(x)} +wm{0,L f V(x)} - \L g V(x)\ (16) 

is negative definite. Observe first that 5(0) = and S(x) < for all x. Assume S(x) = 0. 
Each term of S(x) is nonpositive, so min{0, LfL G V(x)} = min{0, LfV(x)} = \L g V(x)\ = 0. 
By Assumption |2j x = 0, which gives the negative definiteness. Therefore —S(x) is positive 
definite, so we can determine a continuous positive definite real-valued function p such that 
< —S(x) (e.g., p(s) = min{— S(r) : \r\ = s}). We prove that |L g V(o;)| < implies 
that either LfL G V(x) < —\p{\x\) or LfV(x) < —jp(\x\). To this end, consider x such that 
\L g V{x)\ < \p{\x\). Then p(\x\) < — min{0, LfLoV(x)} — min{0, LfV(x)} + by our 

choices of p and S. We deduce that min{0, LfLcV(x)} + min{0, LfV(x)} < ~p(\x\). It 
follows that either min{0, LfL G V(x)} < —\p(\x\) or min{0, LfV(x)} < —\p{\x\). Therefore, 
\L g V(x)\ < \p(\x\) implies LjL G V(x) < —jp(\x\) or LfV(x) < — so we can take 

T(s) = \p{s) and N{s) = \p{s). ■ 
Fourth step. We prove that the right hand side of dT3t is negative definite when the smooth 
positive definite function 5(-) is suitably chosen. By the preceding lemma, there are three cases: 
First Case. \L g V{x)\ < T{\x\) and L f V(x) < -N(\x\). Then the inequality (UB implies that 

U < -lN(\x\) + 5(V(x))L f L G V(x) + £(V(x))8 2 (V(x))\L g L G V(x)\ 2 . (17) 
Choosing 5(-) such that 

S(V(x))L f L G V(x) < ^N(\x\), i{V{x))5 2 {V{x))\L g L G V{x)\ 2 < ^N{\x\) (18) 
for all x G R n . Therefore, CE3> gives U < -\N(\x\) < for all x ^ 0. 

Second Case. \L g V(x)\ < T(\x\) and L f L G V(x) < -iV(|x|). Then the inequalities © and (TT3T) 
imply U < -8{V{x))N{\x\) + £(V(x))5 2 (V(x))\L g L G V(x)\ 2 . Choosing 6(-) such that 

5(V(xm(V(x))\L g L G V(x)\ 2 < 1n(\x\) (19) 

we obtain U < -\5{V(x))N{\x\) < for all x ^ 0. 
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Third Case. |L 9 V(x)| > r(|x|). Then the inequality (IT31) implies that 

U < - l ^(V(x))T 2 (\x\) + 5{V{x))L f L G V{x) + ^V{x))5\V{x))\L g L G V{x)\ 2 . 
Arguing as above provides 5(-) such that 

5(V(x))L f L G V(x) < lT 2 (\x\nV(x)), 6 2 (V(x))\L 9 L G V(x)\ 2 < ir 2 (|*|) , (20) 
so we obtain U < -\l{V{x))T 2 {\x\) < for all x ^ 0. 

Fifth step. To conclude the proof for the control affine case, one has to prove that one can 
determine a C 1 and positive definite function 5(-) simultaneously satisfying the requirements 
(fTTT) . (fl3l) . (EH), (O, (f20l>. This can be done as follows. We can first find a C 1 positive definite 
function 5 satisfying the requirements (fTTT) . (TT81) . (THJ1). d2"0l that is increasing on [0, 1], non- 
increasing on [l,oo) and bounded by 1. We denote this initial choice of S by S a . Next, we 
minorize 1/(1 + A\L G V(x)\) by a positive function of the form x i— > V(V(x)) (using, e.g., 
V(s) = inf{l/(l + A\L G V(x)\) : x E W 1 , V(x) = s}). One can easily determine an everywhere 
positive, non-increasing C 1 function co(-) such that u>(s) < | min {P(s),V(2s), 1} for all s > 0. 
Now consider the function 

Jk, i + a 2 

It is positive definite, of class C 2 , and (since 5 a is bounded by 1) satisfies, for all s > 0, 



\S'(s) 



S a (s)S a (2s)m(s) i 6 a (±s)8 a (s)w(±s) 



, +4? 2 • ^ + Ml*) < • ( 22 ) 

^From this inequality, one can deduce that 5 defined in (|2TT> satisfies (fT3l . On the other hand, 
since to is smaller than 1, the inequality 5(s) < f*, 2 5 a (l)5 a (2l) / (1 + Al 2 )dl is satisfied for all 
s > 0. Now, we distinguish between two cases. First case: If s £ [0,1], then, since 5 a is a 
nonnegative function smaller than 1 and increasing on [0, 1], we get 5(s) < f*, 2 5 a (l)dl < 5 a (s). 
Second case: s > 1, then, since 5 a is a nonnegative function smaller than 1 and nonincreasing 
on [1, +oo), we get 

6 ( s ) < f TZ^2 dl < ^r^ 6 ^ ^ ^ ■ (23) 

Ji s l + s 2 2(1 + s 2 ) 
Hence, the function S defined in (OTT) satisfies the requirements (fTTT) . (fT3l . ( fTSt . dT^l. d2"0"l ). 

Remark 3: The proof of Lemma |2] provides explicit formulae for the functions T and TV 
required for our constructions. On the other hand, the function 5 in © can be obtained by 
simply verifying the requirements in the fifth step of our proof. 
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B. Fully Nonlinear Case 

We now extend the construction to our original fully nonlinear system ©. We can write 

F{x 1 u) = f(x) + g(x)u + h(x,u)u, where h(x,u) = 

Jo 

Along the trajectories of ©, it follows that V = LjV(x) + L g V(x)u + W(x)h(x, u)u. Since F 
is C 2 in u, we can find a continuous function H : [0, oo) x [0, oo) — > (0, oo) that is nondecreasing 
in both variables such that \h(x, u)u\ < H(V(x), \u\)\u\ 2 for all x and u. One can find a A E K,^ 
such that |VV(x)| < a4(|x|) for all x. Taking u to be a feedback of the form © gives 

V < L f V(x) - ^V(x))\L g V{x)\ 2 + H*{V{x), UV{x))L g V{xY\)e{V{x))\L g V{x)\ 2 

with H*(r, s) = a 4 (aj" 1 (r))if(r, s). We now restrict our attention to the set JF[^] of all feedbacks 
© such that £(s) < ^(s) for all s > 0, where we assume the positive function £ is such that 

H m {V{x),Z{V{x))\L g V{x)\)Z{V{x)) < l - VxGM". (25) 

Condition d23l can be satisfied by minorizing £ as necessary without relabelling. (The proof that 
£ can be chosen to satisfy (1231) is similar to the construction of the function 5 in the first part 
of the proof.) Fixing a feedback from this family we get 

V < L f V(x) - l -i{V{x))\L g V{x)\ 2 (26) 

along the closed loop trajectories of (J2j). Applying the construction from the first part of the 
proof to the control affine system (HJ) with £ = £ provides a function 5 and a CLF U of the form 
® such that W(x) := - {L f U(x) - L g U(x)£(y(x))L g V(x) T } is positive definite. Therefore, 
U along the trajectories of © in closed-loop with the feedback ©, with £ satisfying £(s) < £(s) 
for all s > 0, reads U = -W(x) - VU(x)h(x, -£(V(x))L g V(x) T )£(V(x))L g V(x) T . Therefore, 
since H is non-decreasing in its second argument, it follows from our choices of £ and H that 
U < -W(x) + {\VU(x)\H(V(x)^(V(x))\L g V(x)\)£(V(x))}£(V(x)) \L g V(x)\ 2 for all x. One 
can construct a positive nondecreasing function A such that, along the closed loop trajectories, 

U < -W(x) + A(V(x))£(V(x)) \L g V(x)\ 2 . (27) 

Now consider the function © with the above choice of 5 and O(s) = 4A(s), which is positive 
definite and radially unbounded. Then, according to our Assumption [H (l26l) . and (1271) . we get 

< -W(ar) - A(^(x))£(V(x))|L 9 1/(a;)| 2 Vx G R n 



> dF, , 
— (x,An)-— (x,0) 



dA. (24) 
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The right-hand-side of this inequality is negative definite, so we can satisfy the requirements of 
the theorem using fi(s) = 4A(s) and the CLF VK This concludes our proof. 

V. Example 

We illustrate Theorem [T] by applying it to the two-link manipulator (see [1]). This system is 
a fully actuated system described by the Euler-Lagrange equations 

' mr 2 + \ 9 + IMrrQ = r , mr - mr9 2 = F , (28) 

where M is the mass of the arm; L is its length; m is the mass of the gripper; r and 9 denote the 
angle of the link and the position of the gripper, respectively; and r and F are forces acting on the 
system. It is well-known that (l28l) can be stabilized by bounded control laws. On the other hand, 
this system is globally feedback linearizable so a quadratic CLF can be determined. The novelty 
is that we determine a CLF whose derivative along the trajectory is made negative definite by 
an appropriate choice of bounded feedback. Without loss of generality, we take m = M = 1 
and L = y/3. With x\ :— 9, x 2 :— 9,x 3 := r, 24 := r, the system (1281) becomes 

F. (29) 



Xi = x 2 , 



X-2 



2X3X2X4 _|_ T 



x 3 = x 4 , 



X.i 



x 3 xj 



We construct a globally stabilizing feedback, bounded in norm by 1, and an associated CLF for 



d2"9l. We set (p) = 1/(2^1 + p 2 ) for all p e R throughout the sequel. 



Consider the function 

V(x) -- 



(xl + l)x\ + xl + 



1 + x 2 - 2 



(30) 



^3 -r -l )JU2 -r o, 4 -r y 1 + X x 4 

This function is composed of the kinetic energy of the system with additional terms. It is 
positive definite and radially unbounded and its derivative along trajectories of (|2"91 satisfies 
V(x) = x 2 t + X4F + Xi{x 1 )x 2 + x 3 (x 3 )x 4 . Therefore the change of feedback 

r = -xi(xi) + r b , F = -x 3 (x 3 ) + F b (31) 

yields V(x) = x 2 r b + x^F h . On the other hand, after the change of feedback d3lT) . the equations 
of the system take the control affine form x = f(x) + g(x)u with 





x 2 














—2x3x2x4— xi (xi) 




1 







n 


m = 


x^+1 


, g(x) = 


4+1 


, u = 




















_ x\x 3 - x 3 (x 3 ) 




_ 


1 _ 
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Next consider the vector field G{x) = (0, x±, 0, x 3 ) T . Simple calculations yield 

dV dV 

L G V(x) = - — (x)xi + ^ — (x)x 3 = fx, + l)x 2 xi + x A x 3 . (32) 
0x2 0x4 

Since V(L G V(x)) = {x 2 (xl + + l),x 4 + 2xix 2 x 3 , x 3 ), we get 

L f L G V(x) = x|(2x| + l) + xl-a;?(x 1 }-xi(x 3 ). (33) 

We now check that Assumptions [T] and |2] are satisfied. Since LjV(x) = and L g V(x) = [x 2 £4] 
everywhere, Assumption [T] is satisfied. If L g V(x) = 0, then x 2 = x 4 = 0, in which case we get 
L f L G V(x) = -xf(xi) —xl(x 3 ). It follows that if x ^ and L g V(x) = 0, then L f L G V(x) < 0. 
Therefore Assumption El is satisfied. Hence, Theorem [T] applies. Consider the function 

V\x) = A0[2 + 2V(x)} 6 + L G V(x) - 40(2 6 ) . (34) 

Simple multiplications show 80[2 + 2V(x)f > V\x) > 3 (x\ + x\ + x% + x 4 ) for all x, so 
is positive definite and radially unbounded. Moreover, we see that along the trajectories of d29l 
after the change of feedback (|3~T1) . 

V\x) = A80[2 + 2V(x)} 5 (x 2 T b + x 4 F b )+x 2 (2x 2 + l)+x 2 -x 2 (x 1 ) - x 2 (x 3 ) 
+X 1 T b + x 3 F b , 



since V(x) = x 2 r b + x 4 F b . Therefore, from the triangle inequality, we deduce that 

V\x) < s/\ +xjr b 2 + 480[2 + 2V(x)] 5 x 2 r b + x\{2x\ + 1) 



(36) 



+ v /l + x 2 F 2 + 480[2 + 2V(x)fx A F b + x 2 - \x 2 ( Xl ) - \x\(x 3 ) . 
We demonstrate now that V$ is a CLF for d29t by showing that the right hand side of d36t is 
negative definite for the feedbacks 

r b = -x 2 (x 2 ) , F b = -x 4 (x 4 ). (37) 

To this end, notice that we have 

V\x) < Ti{x)x 2 2 {x 2 ) +T 2 (x)x 2 (x i ) - \ [xfixt) + x 2 (x 2 ) + x 2 (x 3 ) + x 2 (x 4 )] (38) 



where we define the TVs by T x {x) = ^1 + x\ - 480(2 + 2V{x)f + 2^/TTxJ(2xl + 1) + \ 
and T 2 (x) = \fl + x 2 3 - 480(2 + 2V(x)f + 2^/1 + x\ + ~. From the expression of V(x), we 
deduce that 7\ and T 2 are nonpositive and therefore 

V\x) < -~ [xKxt) + x\(x 2 ) + x 2 3 (x 3 ) + xl(x 4 )] . (39) 

The right hand side of this inequality is negative definite and the feedbacks resulting from (OH 
and (l3Vb are bounded in norm by 1. This concludes the proof. 

February 2, 2008 DRAFT 



12 



VI. Robustness to Actuator Errors 

Theorem [T] provided a stabilizing feedback u = Ki(x) such that x = f(x) + g(x)Ki(x) is 
globally asymptotically stable (GAS) to x — 0. Moreover, for each e > 0, we can choose K\ to 
satisfy \K x (x)\ < e for all x E W 1 . 

One natural and widely used generalization of the GAS condition is the so-called input- 
to-state stable (ISS) property [18]. For a general nonlinear system x = F(x,d) evolving on 
W 1 x M m (where d represents the disturbance), the ISS property is the requirement that there 
exist (3 E ICC and 7 E K,^ such that the following holds for all measurable essentially bounded 
functions d : [0, 00) — > IR m and corresponding trajectories x(t) for x(t) = F(x(t), d(t)): 

\x(t)\ <P(\x(0)\,t) +T(|d|oo) Vt>0. (ISS) 

Here | • |oo is the essential supremum norm. The ISS property reduces to GAS to for systems 
with no controls, in which case the overshoot term 7(^00) in the ISS decay condition is 0; see 
also [11], [12] for the relationship between the ISS property and asymptotic controllability. It is 
therefore natural to look for a feedback K(x) for (HJ) (which could in principle differ from K\) 
for which 

x = F(x, d) := f(x) + g(x) [K(x) + d] (40) 

is ISS, and for which < e for all x E W\ where e is any prescribed positive constant. 

In other words, we would want an arbitrarily small feedback K that renders (HJ) GAS to x = 
and that has the additional property that (l40l) is also ISS with respect to actuator errors d. 

However, it is clear that this objective cannot be met, since there is no bounded feedback 
K(x) such that the one-dimensional system x = K(x) + d is ISS. On the other hand, if we add 

Assumption 3: A positive nondecreasing smooth function V such that (i) J +o ° ds = +00 
and (ii) |L fl V(ar)| < V(V(x)) for all x E R n is known. 

where V satisfies our continuing Assumptions HE then any feedback K := —£(V(x))L g V(x) T , 
obtained from Theorem [T] for the control affine system x = f(x) + g(x)u and chosen such that 
\^(V(x))L g V(x)\ < e for all x E W 1 , also renders d40l integral-input-to-state stable (USS). For a 
general nonlinear system x = F(x, d) evolving on W 1 x W a , the iISS condition is the following: 
There exist f3 E ICC and a, 7 E /Coo such that for all measurable locally essentially bounded 
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functions d : [0, oo) — > W m and corresponding trajectories x(t) for x(t) = F(x(t),d(t)), 

a(\x(t)\) < /3(\x(0)\,t) + [ 7 (|d(s)|)ds Vt > 0. (iISS) 

The iISS condition reflects the qualitative property of having small overshoots when the distur- 
bances have finite energy. It provides a nonlinear analog of "finite H 2 norm" for linear systems, 
and thus has obvious physical relevance and significance [2], [3]. Assumptions [T]|3] hold for our 
example in the previous section, since in that case, |L 9 V r (x)| < 2(V(x) + 2) for all x G M n , so 
we can take 1>(s) = 2(s + 2). In fact, our assumptions hold for a broader class of Hamiltonian 
systems as well; see Remark |5] below. 

To verify that the Theorem [T] feedback also renders (l40l) iISS, we begin by fixing e > and V 
satisfying our Assumptions [T]|3l and applying our theorem to x = f(x) + g(x)u. This provides 
a CLF U for (HI and a corresponding positive function £ that satisfies \£(V(x))L g V(x)\ < e for 
all x £ IR n . The CLF U has the form ©. By reducing 5 and 5' from Section ITV-Al and replacing 
V(p) with p i — > T>(2p) + 1 in Assumption |3] without relabelling, we can assume 

\L g U(x)\ < V(U(x)) Vx G W n . (41) 

Then 

r u{x) dn 

U(x) = / — f-, where t/(x) = + S(V(x))L G V(x) (42) 

is again a CLF for our dynamic ([l]), since our choice of V gives £/(x) — > +oo as |x| — ■> oo 
because U is radially unbounded, and because VU(x) = VU(x)/T>(U(x)) (which gives the 
CLF decay condition). The smoothness of U follows because U and V are both smooth. Finally, 
ATI) gives 

|L s C/(a;)| = \L g U(x)/V(U(x))\ < 1 Vx G R n . (43) 

We next choose the smooth feedback i^i(x) = — £(V r (x))L 9 V r (x) T , where £ is a smooth positive 
function satisfying the above requirements, so K\ renders (HJ) GAS to x = 0, by Theorem d To 
check that K(x) := K 1 (x) also renders (l40l) iISS, notice that our choice of K x and (1431 give 

Vf/(x)F(x, d) = Vt/(x)[/(x) +^(x)ATi(x)] + L g U(x)d 

< -a 5 (\x\) + \L g U(x)\\d\ < -a 5 (|a;|) + |d| 
for all x and d and some continuous positive definite function a$. Inequality (l44l) says (see [3]) 
that the positive definite radially unbounded smooth function U is an ilSS-CLF for d40b . The 
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fact that (l40b is iISS now follows from the iISS Lyapunov characterization [3, Theorem 1]. We 
conclude as follows: 

Corollary 4: Let the data © satisfy Assumptions [l]|3] for some vector field G : W 1 — > W 1 and 
V : W 1 — ► M, and let e > be given. Then there exist smooth functions 5, £ : [0, oo) — > [0, oo) 
such that (i) the system (O with the feedback := (x))L 9 V(x) T is iISS and has a 

smooth ilSS-CLF of the form <@2|) and (ii) |K(a;)| < £ for all x G R n . 

Remark 5: Assume the Hamiltonian system dU) satisfies the conditions (a)-(b) we introduced 
in Section |TlT| as well as the following additional condition: (c) There exist A, A > such that 
spectrum{M _1 (g)} C [A, A] for all q. (Assumption (c) means there are positive constants c 
and c such that c\p\ 2 < p T M(q)p < c\p\ 2 for all q and p.) Then dU) satisfies our Assumptions 
[l]|3] and so is covered by the preceding corollary. In fact, we saw on p. |5] that (a)-(b) imply 
that Assumptions [l]|2] hold with V — H, and then Assumption |3] follows from (c) because 
\L g V(x)\ 2 = \p T M^(q)\ 2 < X 2 \p\ 2 < (X 2 /X)p T M- 1 (q)p < 2(X 2 /X)V(x) for all x = (q,p). We 
can choose V{s) := ^2{X 2 /X)(s + 1). 

VII. Conclusion 

We showed how to construct control-Lyapunov functions for fully nonlinear systems satisfying 
appropriate generalizations of the Jurdjevic-Quinn conditions. We also constructed feedbacks of 
arbitrarily small norm that render our systems integral-input-to-state stable to actuator errors. 
Our constructions apply to important families of nonlinear systems, and in particular to systems 
described by Euler-Lagrange equations. Redesign and further robustness analysis for our systems 
via our construction of control-Lyapunov functions will be subjects of future work. 
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